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We use a novel method based on the semi-classical analysis of cr-models to describe the phe- 
nomenon of strong localization in quasi one-dimensional conductors, obtaining the density of trans- 
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cases our approach turns out to be exact. The results offer a novel perspective on localization. 
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I. INTRODUCTION 

Beginning with the early work of Thoulessi^, quasi 
one-dimensional conductors have provided a valuable 
arena in which to explore the influence of quantum in- 
terference effects on the transport properties of weakly 
disordered phase coherent conductors. By the early 80's, 
a complete scaling theory of localization in multi-mode 
wires had been formulated by Dorokho^ (and devel- 
oped later by Mello, Pereyra and Kumar 4 ) as a Brow- 
nian motion of eigenvalues of the transmission matrix 
- the 'DMPK equation' (for a recent review see, e.g., 
Ref. H. For several symmetry classes, analytic results 
for the low moments of conductance were obtained in 
both the metallic and strongly localized regimes 5 *^. 
Lately, a third, and potentially more versatile approach 
has been developed to investigate quantum transport 
in (multi-terminal) disordered conductors. By formally 
relating the transmission eigenvalue distribution to a 
'multi-component Green's function' (see below), Nazarov 
has shown£ that known results in the metallic limit can 
be inferred from the equations of motion of the average 
quasi-classical Green function^iiSiii. 

The transport properties of a conductor are fully spec- 
ified by the transmission matrix t 12 . In particular, the 
Landauer-Biittiker formula allows the conductance to be 
expressed through eigenvalues T n of the 'squared' trans- 
mission matrix, 

G =I trttt = & T »- 

n 

The statistics of the T n follow the universal Dorokhov 
distribution, valid in the metallic limit (i.e. where the 
dimcnsionless conductance a = hG/e 2 ^> 1): 

^> e - 2L/£ )-(E^-^)) = f^= T . a) 

n v 

By solving the corresponding quasi-classical equation in 



the diffusive limit, Eq. has been established under 
very general conditions. Yet, while the method intro- 
duced by Nazarov is capable of describing accurately 
the transmission eigenvalue distribution in the metallic 
phase, the method fails to account for the phenomena of 
weak and strong localization. 

In recent years, the development of a field theoretic 
technique to explore phase coherence phenomena in dis- 
ordered conductors has exposed the strengths and limi- 
tations of the quasi-classical scheme. Weakly disordered 
conductors and superconductors can be classified into 
one of ten symmetry classes^, their spectral and trans- 
port properties specified by a field theory of nonlinear 
cr-model typ o 14 ' 15 . It was realized that within the for- 
mally exact framework provided by the field theory, the 
quasi-classical equations of motion represent a saddle- 
point or mean-field approximation, involving only a sin- 
gle field configurationiSiii. Indeed, drawing on the sub- 
stantial literature on the quasi-classical equations in the 
superconducting context (the Eilenberger— and Usadelii 
equations), this identification simplified substantially the 
analysis of physical realizations of the 'novel symmetry' 
classes of disordered superconductors introduced by Alt- 
land and Zirnbauer—. The absence of localization in 
the quasi-classical theory follows then from the neglect 
of fluctuations about this field configuration. 

Describing localization within the cr-model is however 
far from trivial. While weak localization corrections to 
transport can be developed as a systematic perturbation 
theory, the transition to strong localization is signalled by 
the growth of contributions non-perturbative in the con- 
ductance. Indeed, a calculation by Rejaei— of the exact 
transmission eigenvalue distribution relied on a mapping 
of the cr-model for the quasi one-dimensional system onto 
a heat kernel (see also Ref. Il9|) . Such an analysis affords 
little intuition into the origin of localization within the 
framework of the cr-model. In particular, can the transi- 
tion (or, in quasi one-dimension, the crossover) to strong 
localization be understood in terms of certain cr-model 
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field configurations? 

A preliminary answer to this question was provided in 
an insightful work by Muzykantskii and Khmel'nitskii^. 
Using the field theoretic approach, it was demonstrated 
that the long-time response of a weakly disordered quasi 
one-dimensional wire to a voltage step was controlled by 
spatially inhomogeneous saddle-point field configurations 
of the cr-model action. Non-perturbative in the conduc- 
tance, these field configurations were associated with rare 
"anomalously localized" states embedded deep within the 
metallic phase. Whether these states provided a carica- 
ture of states near the Anderson transition remains to 
date the subject of debate. Later, in a related activity, 
it was shown the nucleation of localized quasi-particlc 
states inside the gap of a weakly disordered symmetry 
broken superconductor were ascribed to spatially inho- 
mogeneous saddle-point configurations of the action^. 

In this paper we will show that, for a range of symme- 
try classes, the physics of strong localization is captured 
exactly by the inclusion of a set of new saddle points 
combined with their associated Gaussian fluctuations. 
This recalls a previous investigation by Andreev and Alt- 
shuler— of energy level correlations in symmetry broken 
(zero-dimensional) chaotic systems — random matrix en- 
sembles belonging to the unitary symmetry class. There 
it was shown that the two-point correlator of the density 
of states could be fully recovered from the inclusion of 
two saddle points together with their associated Gaus- 
sian fluctuations. In both cases, the coincidence can be 
traced to the property of semiclassical exactness shared 
by the different cr-model field theories. The aim of the 
present paper is to elucidate this principle in the quasi 
one-dimensional system comparing results for transmis- 
sion eigenvalue distributions to those obtained from heat 
kernel methods. The novel view of localization presented 
here may inform the treatment of more complicated prob- 
lems. 

Following Rejaei, our analysis rests on a relation which 
allows the ensemble average of the generating function of 
the transmission matrix tt^ (at a particular energy) of a 
quasi one-dimensional sample 



I det(l - sin 2 (g/2)ttt) 
\det(l -sinh 2 (0/2)ttt) 



(2) 



to be presented as a partition function of a nonlinear cr- 
model of the appropriate symmetry class. The density 
of transmission eigenvalues may then be extracted using 
the relation 



-i sinh( 



cosh </>„ + cosh ( 



, (3) 



where T n = 1/ cosh 2 (</>„, /2) denote the eigenvalues of the 
matrix ttf . From this function one can infer the trans- 
mission eigenvalue density through the relation 

p{<f>) ee £ (5(4>-4> n )) = l-^+i^-F^-m)) . (4) 



The moments of the transmission matrix can be cast in 
terms of the retarded and advanced Green's functions 

microscopic Hamiltonian 

of the wire as 



tr(ttT)" = tr T {v L G?v R Gf) 

where i>l/r denotes the current operator through left and 
right cross-sections at x = and x = L respectively, and 
the trace on the right hand side runs over spatial coordi- 
nates. Following Nazarov, it is helpful to note that the 
determinants in Eq. can be recast in terms of a multi- 
component Green's function according to the relation 



det(l - 7i7attt) 



tx det 
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det ( * liGfvL 
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where, for example, in the case of the numerator of 
7i72 ee sin (0/2) . Thus the generating function Z(cj>,0) 
may be related to the Green's function of an enlarged 
Hamiltonian, which includes an off-diagonal 'vector po- 
tential'. 

With this representation in hand, the outline of the re- 
mainder of the paper is as follows. In the next section we 
present a complete analysis of the localization properties 
of Class CI, one of the superconducting symmetry classes 
introduced by Altland and Zirnbauer— . We outline the 
cr-model description of this symmetry class and repro- 
duce Nazarov's calculation in this context, before intro- 
ducing the full set of saddle-point trajectories - geodesies 
on the target space - and performing a complete semi- 
classical calculation. From the latter we will extract the 
transmission eigenvalue density and mean conductance. 
The exactness of the results obtained are then verified by 
heat kernel methods. In the following section, the tech- 
niques developed above will be applied to two additional 
symmetry classes that are amenable to such a treatment. 
Finally, in the last section, we will draw conclusions. 



II. CLASS CI 

Superconductors which exhibit both time reversal sym- 
metry and spin rotation invariance belong to symmetry 
class CI. By introducing the familiar Nambu doublet of 
electron and hole operators in what we will refer to as the 
'particle-hole' (ph) space, in the mean-field approxima- 
tion, the quasi-particle Gor'kov or Bogoliubov-de Gcnncs 
Hamiltonian 7i m f can be written as 



ij 
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where and Ay represent real symmetric matrix 
elements corresponding to the non-interacting single- 
particle Hamiltonian of the quasi one-dimensional dis- 
ordered wire and the superconducting order parameter 
respectively. With this structure, the matrix Hamilto- 
nian H satisfies the symmetry relations 



H = H 



-CH'C 



Tn-X 



C 



1 
-1 



(5) 
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which in turn implies the identity 

Gf = -C{G* € ) T C- X , 
between the retarded and advanced Green's functions. 

A. Partition function 

With this definition, the ratio of determinants which 
specifies the generating function J2J) can be conveniently 
expressed via a supersymmetric field integral, with the 
numerator arising from a fermionic integral and the de- 
nominator from its bosonic counterpart. Since these are 
not identical, the corresponding partition function will 
inherit boundary conditions which break the supersym- 
metry. To identify the soft diffusion modes of the super- 
conducting system, it is convenient to effect a doubling of 
the field space of the supervectors ip to include a 'charge- 
conjugation' (cc) spaced 



(40. 



i 



(6) 



Focusing on the transmission eigenvalue distribution at 
e = 0, we notice that if Gq ~ (iptf) is the (1, 1) compo- 
nent of the matrix Green's function Q in the cc space, 
then the advanced component Gq ~ ~C('ip T ?p T )C~ 1 is 
just the (2,2) component. Thus the enlarged structure 



required is just the usual cc space. Altogether, this leads 
to the representation 

det(l — 7i72ttt) 



det(l - Ci<2ttt) 

Z?*X># exp 



g- 1 



where 



n = h <g> i b f ® ice 

+v L T! <g> £+ <g> l ph + v R T 2 <8> E_ 



(7) 
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where Ti.2 = diag(Ci,2j 7i,2)bf > and £j denote the Pauli 
matrices in the cc space. 

Once cast in the form of a functional field integral, it 
is a straightforward (if somewhat lengthy) and standard 
procedure to show that the low-energy properties of the 
partition function are contained with an effective field 
theory of nonlinear cr-model type (for a detailed discus- 
sion of the explicit derivation, we refer to one of the many 
standard references, e.g. Ref »" ,22 ' ). Two points should be 
remarked upon. The first is that, since Class CI has two 
symmetries - particle-hole and time reversal (tr) symme- 
try (see Eq. |5J) - the supervector needs to be doubled 
twice. The size of the supermatrix field Q(r) that ap- 
pears at an intermediate stage in the derivation is thus 
16 x 16 (phxccxtrxbf), but this is quickly reduced to 
the 8x8 field q by Q = a^q on the saddle-point manifold 
(<7j denote the Pauli matrices in ph space), with the ph 
space then disappearing from view. Secondly, following 
Rejaei, the 'vector potential' associated with the coupling 
of the wire to the external leads, can be absorbed into a 
rigid boundary condition on the field integral. Taking the 
metallic contacts at x = 0, L to be to clean normal con- 
ducting leads, the ensemble averaged partition function 
assumes the form (with H = 1) 



J 



det(l 



7i72t e t|) 



det(l - CiC2t £ t!) 



<?(0) = S 3 



X>gexp ( / drSTr[Vg] 2 



(8) 



r 



where the field integral is over 8x8 supermatrix fields Specifically, if we choose 
q subject to the nonlinear constraint <z 2 (r) = 1. Fi- 
nally, the supersymmetry breaking source terms enter 
the boundary condition, through the rotation matrix 



71 = -sin6> , 72 = tan 0/2 



Ci = - sinh ( 



(2 = i tanh</>/2 



S = exp (iT 2 ® £_ ® Itr) exp (iT 1 <g> £+ <g> l tr ) 



then the boundary condition at x = L is q(L) — 
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diag(g b b(£),<7ff(£))bf where 



cosh 4> sinh <f) 
— sinh <$> — cosh < 



COS ft 

isin( 



-zsini 
- cos£ 



it 



(9) 



and we have obtained the required form of Rejaei's rela- 
tion • An alternative approach that matches Nazarov's 
orginal formulation is to use the explicit expression for 
F(<j)) that follows from Eq. © and relate it directly to 
the Green's function in the cc space, i.e. 



F(<f>) 



■ sinh d> ( tr 



tt f (l + sinh 2 (0/2)tt t 



= <tr r [^]>| e=# 

= 5— -4 (trtr {qd x q)H \ x=0 



(10) 



where (• • - ) q denotes an average with respect to the a- 
model action |(HJl, A is the cross-sectional area, and v is 
the single particle density of states. tr tr denotes the trace 
over the tr space. Henceforth, we will work only with the 
one-dimensional field q(x), assuming that the width of 
the wire is much smaller than the localization length. 



B. Nazarov's calculation 

To develop a theory of localization from the field inte- 
gral, it is first useful to establish contact with the theo- 
retical framework introduced by Nazarov. As discussed 
in the introduction, the quasi-classical theory formulated 
by Nazarov is contained within the saddle-point structure 
of the present action: Specifically, a variation of the ac- 
tion with respect to q subject to the nonlinear constraint, 
5 =1, obtains the saddle-point equation 

d x (qd x q) = 0, 

which expresses conservation of the matrix current qd x q. 
Indeed, associating q with the average quasi-classical 
Green function, the saddle-point equation can be inter- 
preted as the generalization of the quasi-classical equa- 
tion derived in Nazarov's early work to the superconduct- 
ing wire. 

The simplest solution satisfying the given boundary 
conditions is block diagonal in the bf space and takes the 
form 



qbb(x) 
qs(x) = 



cosh ^ sinh ^ 
-sinh^ -cosh^ 



1. 



8 X 

COS ^f- 



-1 sin ^f- 



(11) 



If one inserts this solution into Eq. i|lU[l . and neglects 
field fluctuations around the saddle point, one finds 

F{(j>) = -muDA4>/L. 



Then, making use of Eq. Q , one can straightforwardly 
recover the Dorokhov distribution Q for a wire of di- 
mensionless conductance g = 2-kvDA/L. 

However, when derived from the supersymmetric field 
theory, one can see that the saddle-point equation 
presents not just one solution, but a whole family of so- 
lutions in the compact fermionic sector 



cos 
i sin ■ 



(8+2irn) X 
L 

9+2irn) X 



-I SOX 



{6+2-im) X 
(e+2Tvn) X 



n e Z.(12) 



Thus we see that even at the level of a purely semi- 
classical analysis, Nazarov's treatment is missing two el- 
ements; corrections associated with fluctuations around 
the conventional saddle-point solution and so- 

lutions l|12|) that loop multiply around the compact 
fermionic sector. We will show that these two 'channels' 
of corrections are, respectively, responsible for the phe- 
nomenon of weak localization — as, indeed, one might 
have guessed from the connection between the role of 
fluctuations in the cr-model and diagrammatic perturba- 
tion theory — and strong localization. 

To uncover this relation, in the following we will under- 
take the semi-classical analysis taking into account fluc- 
tuations of the matrix fields at quadratic order around 
both the conventional and the non-trivial saddle-point 
field configurations. 



C. Semi-classical calculation 

To perform the field integration over the fluctuations, 
it is necessary to review the structure of the target man- 
ifold. For the minimal (n = 1) nonlinear cr-model asso- 
ciated with the symmetry class CI, the target space - 
spanned by matrices q — wl^u; -1 — turns ou^ to be in 
one-to-one correspondence with G — OSp(2|2), by which 
we mean the Lie group M. + x SU(2) extended to an or- 
thosymplectic Lie supergroup. This group structure of 
the target manifold is the reason for the exactness of the 
semi-classical approximationS&Si. 

The correspondence between q and g £ G is described 
in the Appendix. The space of matrices J!|J) parameter- 
ized by <j) and 9 corresponds to the maximal Abelian sub- 
group A e G 



diag(e*, e" 



l0 ) e A , 



while, in terms of group elements g £ G, the generating 
function is 



g(T)=a{</,,8) 



2(0,6) 



J Vgex V l -± J STr^ 1 g] 2 dt 



ff(0)=l 



(13) 

Here we have switched to the dimensionless variables t = 
x/£, T = L/£, where £ = 2nvDA. (Note that, with this 
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definition, T = 1/g.) With this parametrization, the set 
of saddle-point solutions discussed in the previous section 
corresponds to the geodesic trajectories of free particle 
motion 



= diag(e^/ T , e~^ T , e w ^ T , e~ ie ^ T ), (14) 



where 9 n = 9 + 2im with n G Z. These configurations 
are associated with the classical action 



STr[(a 



W-v-l • (n)-|2 



a. 



r^=(r+^)/4T 



To evaluate the contribution from Gaussian fluctu- 
ations, it is helpful to set g t = Of gt, and express 
the element gt through the exponential parametrization 
g t = cxpXt, with boundary conditions Xq = = Xt- 
Matrices X satisfying the osp(2|2) Lie algebra condition 
X = —eX T e~ 1 with e = dia,g(a x ,ia y )bi may be written 
as 



X 



( d a 0' 

-d 7 6 

S (3 e b 

\— 7 —a c — ey 



(15) 



Then, to get the desired target space with the non- 
compact boson-boson sector Mbb = R+, and compact 
fermion-fermion sector Mg = SU(2), one must take the 
variable d to be real, the variable e to be imaginary, and 
the complex variables b and c to be related by the con- 
dition c = — b. The functional integration measure is 
trivial, T>g = Vg = VX + . . ., up to corrections (of order 
X 3 ) that will not affect the present calculation. 

Using the parametrization for X given in Eq. (|15|) one 
obtains the quadratic action 



The determinants are to be evaluated on the Hilbert 
space of square- integrable functions L 2 ([0,T]) with 
Dirichlet boundary conditions. If z is some complex num- 
ber, the operator D z — ~d 2 + 2(z/T)dt on that Hilbert 
space has eigenf unctions sin(kirt/T) e zt ' T (k 6 N), with 
eigenvalues ((hit) 2 + z 2 ) /T 2 , so its (unregularized) de- 
terminant is given by 

oo 

DetL> 2 = Yl {{^ f + z 2 ) /T 2 . 

fc=i 

Taking the logarithmic differential (to kill the z— 
independent infinity) one obtains 



5\nV)ctD z = 



2z5z 




z * — ' (kir) 2 + z 2 



(—z 1 + coth z) Sz = 5 (In sinh z — In z) 



Therefore, 



„ „ sinhz „ „ 
T)ctD z = x Det_D 



and the above ratio of four determinants gives 
smh(~((f> + i6 n )) sinh {\{<j> — iO n )) fsm9 n 



Finally, when combined with the exponential of the clas- 
sical action and summed over n, one obtains the following 
expression for the partition function: 



S a 



1 



dt STr [X 2 + X[{a { t 



(n)> 



4 dt \d 2 -e 2 -bc- T 



? (be -be) 



-2aS 



T 



'-(ad — aS) 



+2W -^^(7/3 -7/3) 



Performing the Gaussian functional integral over d and 
e is the same as computing the kernel for free particle 
motion on the two-dimensional Euclidean space Lie A, 
and the result is Det _1 (— d 2 ). Performing the remain- 
ing Gaussian functional integrals gives determinants in 
the denominator (bosons: 6, c = —6) and numerator 
(fermions: a, 5 and 7, /3). Collecting all the determinants, 
one finally obtains the total fluctuation contribution to 
the partition function, 



Det 



-d 2 



:,+ i<) 



d t Det 



Det {-d 2 ) Det (~d, 



2iB n 



X ^ sinh( |((ft + i6>n)) sinh(±(0-i6y) 

2^ I, 



■ + iOn) 



h(<t>-ie n ) 



smf 



(16) 



Note that the operator — d 2 + (2i9 n /T)d t starts to ex- 
hibit negative eigenvalues as soon as \9 n \ exceeds ir. This 
means that the Gaussian functional integral over b and 
c = —6 docs not exist in those cases. Thus what we have 
done — evaluating the functional integral in the Gaussian 
approximation around all of the saddle points (without 
paying attention to question of existence) , and summing 
contributions — was a purely formal calculation. Nev- 
ertheless, the answer obtained in this way turns out to 
be exact! Indeed, the validity of this expression can be 
established both directly and indirectly. In the follow- 
ing, we will motivate this conclusion by investigating the 
transmission eigenvalue density and the mean conduc- 
tance from the partition function. Then, by interpreting 
the field integral as a heat kernel, we will find an alter- 
native exact method of computation. 
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D. Eigenvalue density use of the relations ©-0, 

With the partition function in hand, it is a straightfor- 
ward matter to compute the eigenvalue density making 



The first term of this expression is just the Dorokhov 
distribution, and originates from the n = (Nazarov's) 
configuration. In this context, the second term simply 
denotes the weak localization correction associated with 
the n — configuration. All higher-order terms of the 
perturbation expansion in powers of T = 1/q turn out to 
vanish identically, a direct manifestation of semi-classical 
exactness. Indeed, the remaining terms in the result fTTjl 
depend non-analytically on T as e~ 1 / T , and arise from 
the non-trivial saddle-point configurations. 

A non-trivial check on the validity of the expression fol- 
lows from the requirement that pr(0) = for all T > 0. 
(One can infer this condition from the joint probabil- 
ity density given by Brouwer et al&.) The plot of the 
density function pj. l (x) in Fig. ^ shows the characteris- 
tic 'crystallization' of transmission eigenvalues associated 
with the onset of localization. The largest T n (smallest 
4> n ) eigenvalue crystallizes at <p = 2T, corresponding to 
T ~ 4e~ 2T = 4e- 2i /« if L » £. 



E. Mean conductance 

A second, and more direct way of inferring exponen- 
tial localization is to investigate the mean (spin) conduc- 
tance of the superconducting wire. The mean conduc- 
tance (£2 T n ) = C(T) of the thick disordered wire can be 
obtained by integrating T = 1/ cosh 2 (0/2) against the 
eigenvalue density p t ((f>)d(f> with the result 

n—l v 7 

On applying the rescaling C(T) — > 4C(4T), the latter 
agrees with the expression derived in Ref . |y. The leading 
term T~ l is the Ohmic contribution, while the constant 
— 1/3 represents the weak localization correction. 

To extract the large-T asymptotics and demonstrate 
exponential localization, the above expression for C(T) 
is inconvenient. Instead, it is convenient to implement an 
integral transform of the expression which brings it to a 
form in which the large-T asymptotics can be developed. 
Making use of the identities (1+2/a;) e^ 1 ^ = — J (1/x 2 — 
2/x 3 ) e~ 1 / x dx and Xm=i n ~ 2 = 7r2 /6, one may easily 



show that 

/>oo 

C(T) = ]T / ( 1 / t2 - 2^ V/t 3 ) e- n ^ 2 / T dr . 

Noting that the function k i-> (1/t 2 - fc 2 /2r 3 ) e~ fc2/4r 
is the Fourier transform of x i— > 2x 2 e~ x r I^Jttt, Poisson 
summation yields 

Finally, introducing the complement of the error func- 
tion, 

2 f°° _ 2 
erfc(z) = —= / e~ s ds = 1 — erf(z) , 

V 71 " Jz 

the expression for C (t) may be brought to the final form 



C(T) = 




= 4^Zerfc(Zv / T) . (19) 
i=i 

While the original expression provides access to the 
small-T asymptotics, the second one gives easy access 
to the large-T behaviour, viz. 

c(T) T^oo 4 

\JixT 

Recalling T = L/£, the latter shows the characteristic 
exponential dependence of localization with a localization 
length £. 

Consideration of the eigenvalue density and mean con- 
ductance provide strong circumstantial evidence that the 
expression derived for the generating function is exact. 
However, to emphasize the validity of the semi-classical 
expansion, it is useful to provide an explicit calculation of 
the partition function as a solution of the Euclidean-time 
'Schrodinger' equation or heat kernel. 
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4>/T <f/T 

FIG. 1: Transmission eigenvalue densities for the classes CI and Dili at T = 0.01 (weak localization - dashed line) and T = 20 
(strong localization - solid line). Note that although eigenvalue 'crystallization' occurs in both classes, localization in class 
Dili is absent due to the persistence of an eigenvalue near <j> = (or T = 1). 



F. Heat Kernel 

By canonical quantization, the path integral with 
the Lagrangian gSTr (<7 -1 <7) 2 becomes the 'Schrodinger 
equation' for free-particle quantum motion on G. The 
Schrodinger operator is the negative of the Laplace- 
Beltrami operator A (or Laplacian for short) on G. Thus 
the heat kernel satisfies a Schrodinger-like equation which 
is the diffusion (or heat) equation on G: 

d t W(g,g';t) = AW(g,g';t) (t > 0) , 

where A acts on the left argument of W. The reduction 
from the full heat kernel to the radial function u> t (<fi, 9) 
on the maximal Abelian subgroup A takes the diffusion 
equation into 

d t uj t {^9) = (Aw t )(&0) (*>0); 

where we continue to denote the radial part of the Lapla- 
cian by A for simplicity. Since G has superdimension zero 
(there being the same number of bosonic and fermionic 
degrees of freedom) the short-time asymptotics of the 
heat kernel is given by 



+0+ 



(In dimensions d ^ 0, the Gaussian would be preceded 
by a factor (iirt)- d / 2 .) 

To proceed, we must draw on some basic facts from 
the geometry and analysis of the supermanifold G: 

• Diagonalization of matrices g (i.e. g = hah -1 with 
a 6 A) defines a polar decomposition of G. By this 
decomposition, the G-invariant Berezin integration 
measure for G determines a radial integration mea- 
sure Jdxdy, which is positive on the Weyl chamber 
chamber [0, oo] x [0, ir] C LieA 

• By a standard formula™ from the theory of Lie 
groups and symmetric spaces (for a recent pedagog- 
ical review see, for example, Ref. l26j) . the measure 
function J for the case at hand is given by 



J = 



sin 



sinh^ ((4> + id) j2) sinh^ ((</> - i0)/2) 



• Finally, the Laplacian on radial functions is given 

by 

(Aw t )(0,0) = (J-%Jdi + J- x d e Jd B )wt{(l>,6) . 

By introducing complex coordinates z = (4>+i9)/2 and 
z — (<j)—i9)/2 it is apparent that the analytic square root 
of the measure function, 



J 



1 j 2 sin(iz — iz) .cosh(z) .cosh(z) 
sinh(z) sinh(z) sinh(z) sinh(z) 



is harmonic: (9? + dg)J x / 2 — 0. Using this property one 
can easily verify that the radial part of the Laplacian can 
be cast in the form 

(A«t)(0, 9) = J-V^fyfy + dgd e ) J 1/2 uo t {^ 9) . 

It therefore follows that the product E t — J 1 / 2 ^ satisfies 
the Euclidean heat equation 

d t E t = (d\ + d 2 e )E t . 

The Euclidean heat kernel in two dimensions is known 
to be (47rt)- 1 e~^ 2+02 )/ 4t . However, this is not the solu- 
tion we want here: As mentioned earlier, the heat ker- 
nel u>t(4>, 9) is subject to zero-dimensional small-i asymp- 
totics, w t ((f>,0) -> e ~ (02+e2)/4t . This short-time behavior 
is achieved with the choice 



4(9 



K0 2 )/4t 



To confirm that Et satisfies the Euclidean heat equation, 
one uses the identity 9? + <9 2 = d z dz and 



-2Im(z" 



-zz/t 



The solution Q t = J~ 1 / 2 E t thus obtained is not yet 2n- 
periodic in 9 and hence does not lift to a function on the 
Abelian group A = M. + x U(l). Enforcing periodicity by 
summing over images, we obtain: 



E 

nsEZ 



i0 n 



2 + g2)/4t 
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where 9 n = 6 + 2im as before. This is the correct answer, 
and it is easily seen to coincide exactly with the semi- 
classical result for the partition function derived above. 

As emphasized above, the coincidence of the semi- 
classical expansion with the exact expression for the par- 
tition function is a particular feature of the group mani- 
fold structure of the cr-model for the symmetry class CI. 
As a result, we can deduce the existence of other novel 
symmetry classes where the transport properties of the 
quasi one-dimensional system can be inferred from the 
structure of the semi-classical expansion. 



III. GENERALIZATIONS 

Having established the principle, we now turn to con- 
sider extensions of the present scheme to other symmetry 
classes. Here we consider classes Dili and AHI: 



A. Class Dili (even) 

Symmetry Dili is realized^ in superconductors which 
exhibit time-reversal symmetry but where the SU(2) spin 
rotation symmetry is broken by a spin scattering mecha- 
nism such as a spin-orbit interaction. The latter presents 
a critical testing ground for the present theory as it has 
already been established by Brouwer et al. 6 that expo- 
nential localization is absent generically in quasi one- 
dimensional systems of this symmetry class. 

For the symmetry class Dili, the corresponding target 
space of the nonlinear cr-model is a Riemannian sym- 
metric superspace of type C|Di£. This means that we 
are to work again with the complex orthosymplectic Lie 
supergroup, G = OSp(2n|2n), but now one must select 
a non-compact symmetric space Sp(2n, C)/USp(2n) in 



the boson-boson sector and the compact group SO(2n) 
in the fermion-fermion sector. For n = 1 these spaces 
are isomorphic to H 3 (the three-hyperboloid) and U(l) 
respectively. A maximal commuting subgroup A = 
K + x U(l) is still formed by diagonal matrices a = 
diag(e^, e~^, e l \ e~ w ) with e K and 6 G [-7T, tt], so the 
saddle-point configurations (|14fl are unchanged. Only the 
fluctuation contribution differs. 

The parametrization of the target space by the expo- 
nential mapping g = expA is the same as for class CI 
but with the boson-boson and fermion-fermion sectors 
interchanged: 



X 



where d is now imaginary, e real, and c = b. In view of the 
duality (by exchange of the compact and non-compact 
sectors) connecting the nonlinear cr-models for the classes 
CI and Dili, all calculations for Dili are very similar to 
those for CI and, for brevity, we simply quote the results 
here. The partition function, obtained by performing 
the sum over geodesies together with the integral over 
Gaussian fluctuations, is given by 




sinh(§(0 + i0„)) sinh (|(</> - i6 n )) 



n6Z 



.) 



sinh ( 



.-O +0/ 4T 



(20) 



From Rejaei's relation, the corresponding density of 
transmission eigenvalues is then given by 



Dili 

Pt 



(0) 



1 

2T 



+ 7T 2 ) 



n#0 



-n(n+l)ir 2 /T 



2lT 2 7 



Re 



cf> + in(n + 1) 



r 



(21) 



A plot of this function in Fig.^shows the 'crystallization' 
of transmission eigenvalues for t ^> 1. Yet, exponential 
localization does not take place, as Pt{<I>) peaks at 4> — 
(maximal transmission), and the peak amplitude decays 
only algebraically with increasing T . 

Turning to the mean conductance, the partition func- 
tion yields the following expression: 



C(T) = 



T 



E 

n=l 



-n 2 n 2 /T 



The first expression agrees with the result of Ref. |(| the 
second is obtained from it by Poisson resummation. The 
first two terms, C pcrt (T) = t" 1 + 1/3, represent the 
Ohmic and weak anti-localization terms that can be com- 
puted from a standard perturbation theory for the mean 
conductance. As with symmetry class CI, higher-order 
corrections from the perturbative expansion vanish iden- 
tically. All non-zero corrections are non-perturbative and 
arise from the non-trivial geodesies in the compact sector. 
The second expression in Eq. (1221 presents an anomalous 
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diffusive asymptotics for large T = L/£: 




(T » 1) , 



a feature already seen in Ref. 0- 

Finally, expressed in the form of the heat kernel, the 
solution for class Dili is trivially related by the duality 
discussed earlier and the validity of the expression for the 
partition function above may be confirmed straightfor- 
wardly. As a final application of semi-classical exactness, 
we turn now to one of the chiral symmetry classes. 

B. Class AIII 

This symmetry class is relevant to the low-energy 
physics of the chiral Dirac operator— and to the ran- 
dom flux model 28 . Lately, it has been discusse d 2 ?^ 
in the context of the quasi-particle properties of a d- 
wave superconductor subject to a smooth random poten- 
tial. Previous studies of quasi one-dimensional systems 
in this class2ii22i2&2i have revealed a surprisingly rich be- 
haviour, including localization-delocalization transitions 
with a number of critical points. 



where V m projects onto the 2m site of each chain. It 
seems that these moments cannot in general be presented 
as the coefficients in some expansion of a functional deter- 
minant. The average conductance (n — 1) may however 
be obtained as the <Pl<Pr coefficient of the determinant 
of a hamiltonian with hopping e"^ i between the 1 st and 
2 nd sites, and between the 2M - 1 th and 2M th sites 
(since hopping is described by a Sj term in the Hamil- 
tonian) . 

We may thus study the usual partition function 

g T =a(<p,e) / T \ 

Z T (<P,9)= f Vgexpl— / ' Cdt\ , 



for a((f>,9) — diag (e l< ^, e 10 ) (the matrix g being 2 x 2). 
The conductance can be then identified with the coeffi- 
cient C(T) in the expansion 

Zr(<t>,0) = l + C {T){cp-ie)~C{T){<jP + 9 2 )/A 
+C 1 (T)(<b-i6) 2 + ... (24) 



1. Rejaei relation 

Before exploring the properties of the nonlinear a- 
model action, it is first necessary to confirm the form of 
Rejaei's relation for this symmetry class. For this we refer 
to the cr-model formulation of the bond disordered quasi 
one-dimensional chain given by Altland and MerkljM. We 
consider N chains of 2M ^> N sites with hopping ma- 
trix element = 1 + (— ) m a along the chains between 
sites 2m and 2m + 1, where a is a 'staggering' parame- 
ter. In the following, we suppose that a weak disorder 
potential (which, further, lifts time-reversal symmetry) 
of strength -C 1 couples the chains. We view the single 
particle Hilbert space as the product H = C 2 ®C N <g>C M , 
with the first factor referring to the odd and even num- 
bered sites along the chain. 

The sublattice symmetry of the Hamiltonian translates 
to the condition H = — H^H^^, where Ej represent the 
Pauli matrices in this odd/even subspace. From this sub- 
lattice symmetry there follows the relation 

G\ = -E 3 G? E 3 . 

The corresponding current operator is simply given by 
E2/2 (if we ignore the disorder and staggering, or evalu- 
ate it in the clean, unstaggered leads). Thus the moments 
of the transmission matrix at zero energy are given by 

tr(totj)" = tr H (v L G$v R G§) n 

= tr H (PiE 2 G^ M E 2 G^) n 

= — trn {V\ Ei Gq'Vm Ei Gj^) n , (23) 



2. Semi-classical calculation 

The most general nonlinear cr-model Lagrangian for 
quasi one-dimensional systems belonging to class AUI 
(and in suitable units of length) depends on two param- 
eters u, v. 

C = STr{g- 1 gf + AvSTr g^g - uiSTrg^g) 2 . 

The target space of the nonlinear cr-model for this classic 
is of type A\A and is obtained from the complex Lie su- 
pergroup GL(n|n) by picking Mbb = GL(n, C)/U(n) in 
the boson-boson sector and Mg — U(n) in the fermion- 
fermion sector. For n = 1 these are GL(1, C)/U(l) ~ K+ 
and U(l) respectively. The exponential parametrization 
g = exp X of the target space is achieved by setting 




with real commuting variables p, q and anti-commuting 
variables a, 0. Again, a(<f>,6) parameterizes the maxi- 
mal Abelian subgroup. Based on the geodesies a[ n ^ = 
diag (e^ t/T , e i9nt/T ) with 9 n = 8 + 2nn, the semiclassical 
approximation for Zt(4>,0) is found to be 
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= E(-D r 



sinh(§(^-i0)) 



n£ 



-v^—iQ^ 



2 +el-u(<p-ie n ) 2 ]/4T 



(25) 



r 



Expanding the result above for the partition function, 
and using Eq. gives 



C(T) 



T T 



£(-!)" cos(27 



-(l+u)n 2 7r 2 /T 



y/nT(l + U) 



^ e -(l+v-l/2) 2 T/(l+u) ^ 



where the second expression is obtained by Poisson re- 
summation of the first. 

The result (|26l) depends periodically on the parameter 
9 = 2?tv with period 2tt. In fact, has the meaning of a 
topological angle. From Ref. |3J we know that tuning 6 
from to 2-7T amounts to changing by two the number of 
channels (which is assumed to be large in order for the 
a- model approximation be valid) . In terms of the bond 
disordered chains discussed earlier 

9 = n(N-f) mod27r, 

where / is the staggering parameter a scaled by the 
strength of the disorder. We see that without staggering, 
a delocalized state appears only for odd channel num- 
ber. The localization length diverges with critical expo- 
nent 1 as the staggering approaches the critical values 
f = N — 1 (mod 2). 

If class Alii is realized by Dirac fermions in a random 
Abelian gauge field, the parameter u > measures the 
strength of the gauge disorder. Although the Lagrangian 
C at first sight would seem to become unstable at it = 1 
(the coefficient of p 2 becomes negative there) , this is not 
really so. What matters is the full quadratic form 



r 



(1 — u)p 2 + (1 + u)q 2 + 2iupq + odd variables 



which has the Jordan normal matrix form 



1 u 
1 



and is 



invertible for any u. In fact, the final result l|2t)l) reveals 
a smooth dependence of the mean conductance on u for 
all physical values (u > 0) including the fake singularity 
U = 1. 



3. Heat Kernel 

Finally, seeking an exact solution through the heat ker- 
nel, an exact calculation of ujt{4>> &) is still possible as the 
Hamiltonian TL obtained by canonical quantization from 
C can once again be transformed to Euclidean form: 

H = -J- 1/2 (fyfy + d e do + 2v(d 4> - id e ) 

+u(d^ - ide) 2 ) J 1/2 , 



with J = 1/ sinh 2 {\{4> — i0)) . By solving the differential 
equation dtuJt — ~Hu>t with the appropriate (5-function 
initial conditions at t = 0+, we exactly recover the semi- 
classical answer (I25|l . 



IV. CONCLUSION 

We have confirmed the exactness of our semiclassical 
analysis for symmetry classes CI, Dili, and AIII. By sim- 
ply summing over the saddle points of the classical action 
— recall that these differed by the length of the geodesic 
looping around the compact fermion-fermion submani- 
fold of the theory — and treating the fluctuations in the 
Gaussian approximation, we have obtained the exact re- 
sult for the partition function. 

The underlying reason for the exactness of our calcula- 
tions is that the path integral on the cr-model manifolds 
considered here satisfies an infinite-dimensional general- 
isation of the Duistermaat-Heckman theoremSiSi. That 
is, the integration manifold in Eg. 1131 being related to the 
space QG of based loops in G, is symplectic with Liou- 
ville measure, and the action is the momentum mapping 
on this manifold. 

We have discussed all the symmetry classes where the 
cr-model target space has a group-like structure, and thus 
exhausted the situations in which the semiclassical ap- 
proach is exact. Nevertheless, we hope that this approach 
may inform future investigations of more taxing localiza- 
tion problems. One very interesting extension that could 
be tackled in a similar manner is the critical scaling of 
the localization length near the band center in situations 
where zero energy is described by class Z3III or ALU. This 
lies outside of the scope of methods like the DMPK equa- 
tion. 



APPENDIX: EMBEDDING OF g IN Q FOR 
CLASS CI 



Let us review the structure of the cr-model targ et space 
for Class CI. The construction used in Ref. [l3| is based 
on 'copying' the symmetries of the Hamiltonian to sym- 
metries of the auxiliary space of the Q- field. For Class 
CI (particle- hole and time reversal symmetry), this is 
achieved by demanding that the fields satisfy 



*M = 



(A.l) 



The subscript M is to remind us that in the formulation 
of Ref. 1 13 the variables and are (super) matrices 
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mapping the auxiliary space to the physical Hilbert space 
(where the Hamiltonian acts), and vice versa. In the 
present case we work with the more traditional column 
and row supervectors VfV and ^v, since: 

• In the localization problem the Q(r) field used 
to decouple the term arising from potential dis- 
order retains the Hilbert space structure of the 
Hamiltonian (particle hole space and spatial index), 
whereas in the random matrix problem that struc- 
ture is lost on averaging and, more importantly, 

• We use some of the additional structure (the cc 
space) to enlarge our Hamiltonian (Eq. 01 • 

Transcribing the conditions l|A.l|) for supervectors gives 



(A.2) 



where we used the fact that 7 and r can be chosen to 
be orthogonal matrices. Now in the localization problem 
Q(r) is a 16 x 16 matrix (phxccxtrxbf - by 'tr' we mean 
the space to accommodate time reversal symmetry), so 
we need to check that the conditions (|A.2|I give rise to a 
cr-model with the same 8x8 space as in the random ma- 
trix problem. This is straightforwardly verified following 
Ref. by observing that Q(r) ~ $ ® ^/a 3 (we drop the 
V subscript from now on), arising from a decoupling of 
potential disorder, has the symmetry 



Q = a- 3 TQ T T~ 1 a 3 . 
The saddle-point manifold is spanned by 

Q = Wa 3 j: 3 W-\ W = w®t ph , 



or 



Q = <73<?, 



q = wL, 3 w 



Thus (|A.3(1 implies that the symmetries of the 8x8 field 
q are 



T -1 

q = -jq 7 , 

T -1 

q = rq t , 



(A.4) 

which are the same relations found in Ref. 0. To con- 
struct an embedding of the 4x4 matrix group-valued 
field g in q requires an explicit choice of 7 and r. We use 

7 = ^bb ® 7b + Eg <X> 7f, 
T = E hh ®T b + E S ®Tf, 



7b = 
n, 



Si 

= a. 



r 3 , 7f = 1^2 6 

<8>Tl, T{ = S3 



Itr, 
5 IT 2 , 



where Tj are the Pauli matrices in tr space, and £"bb and 
Eg are projectors in the bose-bose and fermi-fermi sec- 
tors. Note that the structure of the vector in cc space in 



the bosonic sector is in fact then not given by Eq. H3 but 
rather looks like 



Fortunately this doesn't change the boundary conditions. 
Now the correct embedding is obtained by finding a 
transformation to diagonalize the matrix 

77 = -i7T _1 = El (g) T2 ® Ibf, 



U*T]U = 



(1 

10 

0-10 

Vo -1/ 



® Ibf, 



while simultaneously sending 



Itr ® Ibf 



0-10 
-1 
-10 
0-100 



<g> i b f. 



Why? The point is that the set of w's that preserve the 
symmetries (|A.4|) takes the form 



diag(<?i,<? 2 ), gi ,g 2 e OSp(2|2), 



(A.5) 



in such a basis. The stability group of transformations 

(A.3) 

that send wT, 3 w 1 — * £3 is then given by Eq. I A. 51 with 
.9i = .92, so that the saddle-point manifold is parameter- 
ized by 

q = C/diag( 5 , g- 1 )^^, g S OS P (2|2) . 
An explicit form for U is 



i -1 i 1 

i 



U 2\-i I 1 



® ibf- 



1 



1 



Let us verify that this takes the maximal abelian sub- 
group A to the q matrices of the form 0: 



q = Udi&g{a{4> 1 6) 1 a- 1 {(t) i e))U^ 3 

I cosh 6 — sinh cp\ ( cos f 



diag 



ysinh ( 



cosh < 



i sin 8 — cos 9 



bf 



Finally, the form of the action i|13|) in terms of g follows 
from 

STr[Vg] 2 = 2STr [VgVg- 1 ] 

= -2STr[«?- 1 V 5 ] 2 . 
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